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SUMMARY 


This contract mid-term report discusses the progress of development of 
analytic procedures and computer codes for prediction of weight optimized 
radioisotope thermal' generator shields for unmanned spacecraft operating 

’in vacuo’. Optimization code S0SC, designed to determine shield 

optimum weights and dimensions with respect to specified criterion fast 
neutron plus gamma photon fluxes, is described. The code employs a 
combination of analytic, albedo and Monte Carlo techniques. A theoretical 
discussion and example predicted shield data are given. 

The work effort required to complete the contract program is reviewed . 



1 . INTRODUCTION 


This mid-term report, prepared for the National Aeronautics and Space Ad- 
ministration, Goddard Space Flight Center, by NUS CORPORATION, under 
contract NAS5 -11649 , discusses the progress of development of analytic 
procedures and computer codes for prediction of weight optimized radiation 
shields for an unmanned spacecraft operating 'in vacuo' . 

The work effort has been specifically oriented to the radiation field at the 
energetic particle experiment package indicated in the spacecraft configur- 
ation of Figure 1 . Spacecraft general dimensions , deployment distances 
and materials were obtained from preliminary design drawings furnished by 
NASA-GSFC. The work has considered the radiosisotope generators (RTG's) 
as being plutonium -oxide fuelled, viz. the SNAP-27 RTG. Although the source 
gamma photon and fast neutron angular-energy spectrum will be furnished by 
NASA-GSFC, a number and energy distribution has been assumed based on 
SNAP-27 and PuC >2 published data. 

A shield optimization study code S0SC was designed to determine the 

material thickness and weight required to limit the spacecraft mission ex- 
periment package to a radiation flux exposure of ^ 10 particles/cm^-sec . 

The incident flux was considered as being the -sum of gamma photons and 
neutrons either transmitted by the shield or scattered by the spacecraft 
structure . 

Code S0SC will predict shield requirements for the case of gamma photons, 
according to a combination of analytic particle transmission theory, the 
Monte Carlo transport .method and the albedo technique (backscattering 
theory). It will employ three component sub-codes for this determination. 
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namely: XEST , NUGAM1 and ALB. Preliminary versions of these codes are 
presently designed, debugged and operational. 

Although code S0SC is provisionally designed to evaluate fast neutron trans- 
port in a manner similar to that for gamma photons, it was concurred early 
in the present contract work to focus on photon transport in this stage of 
the NASA program. The code is presently designed to evaluate neutron 
transport using relaxation theory methods. This course for the case of the 
SNAP-27 is substantiated by the fact that the RTG total neutron emission 

n 

rate in the axially perpendicular direction is reported as being 5.7x10 
n/sec^ . This is in good agreement with, but less than an earlier NUS 

Q 

estimate of 1.0x10 n/sec (progress report #1 attachment) . Taking the 

gamma photon dose rate as one-tenth of the neutron dose rate'and allowing 

for dose-to-flux conversion as well as spectral distribution gives an inte- 

9 

grated RTG emitted photon source strength of ~1.0xl0 y/sec, or approxi- 
mately ten times that for neutrons. 

Section 2 of this report describes the theory underlying code S0SC, a 
description of the code and some observations resulting from preliminary 
data determined by the code. Section '3 discusses the work program for 
the remainder of the contract period. 
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2. CODE DESCRIPTION 


2 . 1 Introduction 

This report section describes the general theory, logic and discusses pre- 
liminary results obtained using S0SC , a code written in the FORTRAN-IV 
language for the NASA-GSFC IBM-360/91 digital computer. Code S0SC 
presently consists of three provisionally operational component codes: 

XEST, NUGAM1 and ALB. The reference here to provisional means that 
although these sub-codes are presently operational, they are not yet tested 
to their complete designed capability. Code XEST predicts approximate shield 
dimensions based on analytic methods. Code NUGAM1 predicts shield build- 
up factors for final shield dimensions, by the Monte Carlo method. Code 
ALB determines flux intensities resulting from scattering by the spacecraft 
structure, by means of the albedo technique. 

2 . 2 Theory 

For the purposes of this section the complex spacecraft configuration shown 
in Figure 1 is redrawn schematically in Figure 2. Only one RTG source is 
indicated and the mission experiment package is referred to as a detector. 

In addition, the spacecraft body is replaced by a simpler geometry and the 
boom-arms omitted. 

The radiation number flux at a detector distant r Q from a source S (E 0 ), of 
neutrons or gamma photons E Q , as in Figure 2 , without a shadow shield , 
may be defined as 

0D (E d ) = 0a (Ea) + 0 O (E Q ) , (1) 


where 



< pa (Eo^ - 0ois (Eos) + 0o;R (Eo:r) 


<pa.s (Eois) = the primary radiation number scattered to the 
detector by the area A composed of material j , 
distant r^ and r 2 from the source and detector, 
respectively, 

0aR (EofO = the 'area A ‘-originating reaction product flux 
reaching the detector , 

= 0, for gamma photons as the primary radiation, 
excepting photoneutron interactions (y,n ) , 

0o (E 0 ) = the uncollided primary radiation number reach- 

ing the detector. 

The energy arguments signify that the detected flux consists of radiation 
of primary energy E . scattered energies Ea s (<£<-,) and reaction product 
energies, E R . The subscript a refers to the spacecraft structure as a secondary 
source, eg. scattering, i.e. to albedo fluxes. 

The number flux reaching the detector with a shadow shield,composed of 
material i, as indicated in Figure 2, may be defined as 

0D (E d ) = <pa (Ea) + 0^ (E Q ; E a ) , (2) 

where 

0^ (Eo; E a ) = the shield attenuated flux, 

= (E 0 ) + 0 3 (Eg) , 

0k (E 0 ) = the number flux transmitted by the shield 

without interaction. 
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$a ^a^ 9^s ^ E as^ + ( E ap) + 9^R ( E aR^ < 


^s (^as) “ the shield forward scattered radiation number 
flux reaching the detector, 

<^kp(E a p) = the gamma photon number flux resulting from 
pair production interactions in the shield, 

= 0, for incident photon energies E c ^1.02 MeV, 

= 0 , for incident neutrons , 

(E 3 r) = the shield -originating reaction product number. 

flux reaching the detecto^ e.g. gamma photons 
resulting from (n,^) interactions, 

= 0, generally, for gamma photons as the primary 
radiation, excepting such as photoneutron inter- 
actions , i.e. (y,n) . 

For fast neutrons as the primary incident radiation, the number flux terms 
0 o;r (E^r) anc ^ ( E aR^ Equation (1) and (2) refers to all product radiations, 
eg. neutrons, gamma photons, alphas, protons, depending on the reaction 
probabilities for each. 

The number flux terms in Equations (1) and (2) may be estimated either from 
a combination of analytic relationships and published empirical data or from 
experiment, either numerical analogue, ie. Monte Carlo method or the con- 
ventional laboratory kind. The sole use of the Monte Carlo method is con- 
sidered as being uneconomical and unjustified. A laboratory experiment 
is planned for the future by NASA-GSFC as part of the overall program. 

The present work is thus confined to number flux predictions obtained by 
analytic methods and by published empirical data and judicial use of Monte 
Carlo techniques . 
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Neglecting the reaction product number flux terms for the present, the 
detector incident number flux may be written as 

0D (Epj) = ( E o) + 0b ( E a ) + 00's ( E as) . (3) 

For a normally incident parallel radiation flux 0£,(E O ) , the uncollided number flux 
transmitted through a’ shield of thickness L and reaching the detector, is obtained 
as 

■fe <>y - 0 O < E o 5 e ~ M 0 (4 > 

where 

H (E 0 ) = the total linear attenuation coefficient of the shield 

material for radiation of energy E Q ; the 
notation £ is generally used for neutrons. 


The ratio of the shield total -to -uncollided transmitted flux is referred to as 
"build-up" . The build-up factor may thus be defined as 


R t } _ + 0fo (5 ;_i) 0r, (EntE;.,) 

M o,L) 0a (E 0 ) 0a (Eo) • 


(5) 


The total shield transmitted flux at the detector may be obtained from 
Equations (4) and (5) as 


^ <E 0 ;E a )= B (E 0 , L) . 0 O (E c ) . e‘F E °>- L , 


( 6 ) 


if the build-up factor is known. 

(2 3) 

Although energy and dose build-up factors may be obtained for gamma photons' ' 
and to a lesser degree for fast-neutrons they pertain in almost all cases to semi- 
infinite single -material-composition shields. For small finite shields, i.e. shadow 
shields, and number flux requirements as opposed to energy and dose, recourse 
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to either a laboratory experiment or a Monte Carlo study is a prerequisite . 

In the present work a Monte Carlo evaluation is underway for gamma photons 
and proposed for fast-neutrons. As noted in Section 2.4 , the number build- 
up factor for photons was predicted as: 1 sB s 2, for )|L< 4, even for low 
atomic number materials . 


For a spectrum of incident source particle energies the shield transmitted 
flux is obtained by integration as 

l E a> -/” &< E o- E a> dE ° W 

q 

“2 0a ( E ok' E a ) A E k (8) 

k=l 


For a stratified or homogenous shield composed of m materials, each of thick- 
ness 1 Equation (8) may be rewritten 


q 

$ a (E a ) - S 
k=l 

m 

0o ( E ok) JJ s ( E o ’ 

m 

£.) e “Mi ( E ok)^i 

(9) 

q 

= s 

k=l 

cfio (Eok). B 0E o ) m 

e ~W. ( E ok) *^i 

(10) 


where 


B (E Q ) m = the build-up factor for the composite shield 
of m materials and incident radiation of 
energy £ Q ., and given geometry. 


i = material identity index. 


The uncollided number flux in the foregoing equations may be determined 
from a relationship of the kind 


0 o ( E o> - 0o ( E o' r o } = S < E o> * G < r o> 


( 11 ) 
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where 


S (E ) = source emission rate for radiation of 

energy E Q , 

G (r Q ) = the geometry relationship for a source-to- 
detector distance r Q . 

For example , for an isotropically emitting point source and r Q » detector 
lateral extent, the geometry factor is 

G <r Q ) = (4 ir r o 2 ) _1 (12) 

Extended source and detector geometries may be evaluated according to 
the 'Point-Kernel Method 2 ^ . 

The number flux term , 0 aS (E^g) in Equation (3) , resulting from primary gamma 
photons scattered by an area A, as in Figure 2 , may be redefined as 

= Ja 0 as (E o' Q ' 6 o' <pt V r 2 7 t; ^ dA/ (l3) 

where' 

dA = the differential scattering area, 

0 O = the angle between the incident radiation direction 
and the outward normal of area dA , 

0 = the angle between the emergent (scattered) radia- 

tion direction and the outward ncrmal of area dA, 

{ p - the azimuth angle of scattering in the plane of 

area A, 

r-^ = the distance between the source and the area dA, 

rg = the distance between the area dA and the detector , 
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t = the thickness of the scattering material at area dA , 
measured along the inward normal to dA, 


i = the identity index of the material of which dA 
is a part. 


For a spectrum of source particle energies the energy integrated flux, 

# c (E _) ., may be obtained by an integration similar to that of Equations 

0 !® CXP 

(7) and (8) , as 

00 

V <v =£ <V (14) 

The flux term in the integrand of Equation (13) may be defined according to 
albedo theory as^ 2,4 ^ 


0 aS ( E aS ) = 0 o (E o ,r 1 ) . , cos ?o • « a (E 0 ,9 0 ,Q,(p,t;i) , 


r 2' 


(15) 


where 


a(E 0 ‘Boid ,cp,t;i) = the angular differential number current 

albedo with respect to the noted argu- 
ments (defined for Equation (13))^ 2 ' 4 ^ , 

0 O (E q ,r^) = number flux incident on area dA, 

= S (E Q ) . G (r-^ , c.f. Equation (11). 

The assumption underlying the use of the albedo technique for complex geo- 
metry analysis is that the scattered radiation particles emerge from the 

scattering medium surface at a point close to their point of entry. This 

(2 ) 

assumption is generally justified , eg. the separation distance between 
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entry and exit for one-half of all escaping gamma photons has been found 
to be less than one mean-free-path (for incident energy E 0 ) . Photon 
scattering from very thin or laterally small structures of volume V, may be 
alternately predicted by the single -scattering approximation method^ 5 ^ , 
from the relationship 


0(ySS / 


( E o< 
Volume , 


V 


■l) 


N c 


ctkN (Eq. 



9 a ) dV 


(16) 


where 


N e = the scattering material electron density per 

cubic centimeter, 

a KN^ E o'®s^ = Klein-Nishina angular-energy intensity 

( 2 ) 

distribution function , for photons of 
energy E Q and scattering angle 0 S , 

0 S = the angle between the primary and scattered 

photon directions . 

Equation (15) may be solved' if values for the number albedo are known. The 

albedo may be determined by either a laboratory experiment or a Monte Carlo 

treatment . For gamma photons , recently developed modification of the moments 

method has been reported as a potential source of albedo data' ' . In the 

(7 8) 

present work, experimental albedo data' ' , and empirical relationships 

(2 4 9 10) 

in accord with both experimental and Monte Carlo results • • • , were 

used for gamma photons. A similar approach is proposed for the case of 
fa st -neutrons . 

Since weight is the product of volume and density , the weight optimization 
of an axially symmetric shadow shield of specified composition may be con- 
sidered as an optimization of shield thickness, L m £ n , such that 
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9 

(E^) y , n <sC, where C is a specified criterion, eg. 10 particles/cm -sec. 
The optimum , or minimum weight of a right-cylindrical shield of radius R , 
may be obtained as 


^min 

2 

= ttR L m i n . p , 

(17) 


2 m ' 

- irR S Pi » 

i=l 

(18) 


where 

= the density of shield material i, 
p = the weighted density of the shield / 

= the shield material actual density if m = 1 . 


The total number flux reaching the detector for the case of a polyenergetic 
source, is obtained from Equation (2) doubly integrated over primary and 
secondary energies, as . 




Ld 




~Je q ,E 0o; (E ^y'n/ dE o dE O! + jE o ,E^ (E o ;E a ) y dE o dE o! + 


+f E ( E 0'" E a) n dE 0 dE^ (19) 

where the subscripts y and n denote gamma photons and fast neutrons. The 

use of the albedo and build-up factor concepts is tantamount to an integration 

over E and E , respectively. 

• a a ■ - 


Equation (19) may be further redefined as 



( E o , ' E a)‘y dE o dE a + Je E 9^( E 0' E a) n dE o dE a 1 


( 20 ) 
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where 


(E dV' n^D'/Eo^^V y 'n dE o dE ft 
= C ~Je o/ E^ E o:V' n dE o d E a 


( 21 ) 


( 22 ) 


The flux terms in Equation (20) may be obtained from Equations (10) and 
(13) and L m in obtained by an iterative solution. For example, for a. single 
material shield exposed to a monoenergetic source of neutrons and photons , 
Equation (20) reduces to 


p - 0O < E 0> r B <E 0 ) y e'M®0 ^ (E Q ) n B (E 0 ) n • L , ( 23 ) 


The foregoing theoretical discussion has presumed a knowledge of gamma 
photon and fast neutron interaction phenomena. Such phenomena and the 
relevant interaction physics are summarily reviewed in the Appendix. In 
addition, a familiarity with the solution of radiation transport problems 
by means of the Monte Carlo method is presumed; the reader is referred 
to the references in this regard^ 1 ’*''. 12 ^ . 

The second term in the right-side of Equation (-20) is' defined by Equation (2) . 

It includes the shield originating neutron reaction product flux <^&R(E a R) , 
as yet not discussed in any detail. Although the flux 9 ^r(E 3 ^) reaching the 
detector, may be predicted by means of either a laboratory experiment or 
a Monte Carlo code evaluation, it may be estimated for the case of reaction 
product gamma photons such as result from inelastic scatters or absorptions 
in a'n axially symmetric shadow shield, as 

■ &R< E aR> =/ 0 L g.G(r t ).^(Ep.^(Eo) n 41-e' XT ).e _a ° £(E ° U . e '' i( ^R )(L “ <) dt, (24) 
where 


g = the cross-section area of the axially symmetric 

shadow shield. 
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I 

G(r t ) 

r t 

^ot (E o^ 

0o^ E o^n 

|i( E aR) 

X 


= ttR^ , for a cylindrical shield , 

= the distance from the shield face at the source 
to the differential volume g.d-t, 

= the geometry factor for distance r c.f. 
Equation (11) , 

= the distance from the differential volume 
g . d t to the detector , 

= the linear attenuation coefficient of the shield 
material for fast-neutrons of energy E Q , for 
production of gamma photons , 

= the total- linear attenuation coefficient of 
the shield material for fast-neutrons of 
energy E Q , 

= the shield normally and parallel incident flux 
of fast neutrons of energy E Q , 

= 0 o ( E o'( r o" r ^)n' c.f. -Equation (11) , 

= the total linear attenuation coefficient of the 
shield material for reaction product gamma 
photons of energy E a ^ ; 

= the radioactive decay constant of the reaction- 
produced or compound nucleus , 

= the duration of exposure to the neutron flux , 
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For inelastic scattering the decay constant in Equation (24). is relatively 
large and thus 

1 - e _XT - 1.0 , - (25) 

This is also true for activation where the product \.r is large. 

A modified form of Equation (24) , where g ,Z,r^ and the exponent Ju(E s r) (L~l) 
are replaced by dA, t, ^ and p (E^r) respectively, may be defined as 

2 >..S/E 0 ) .0 o (E o ,r 1 ) n (l-e' XT )e' 25 r°t (E ° )t .e'» l(E *R )t dt (26) 

J ° 

to estimate the spacecraft structure reaction product 'albedo' , for 
prediction of fluxes , 0 o ,r(E q r) . This albedo , valid for normal incidence 
and emergence, may be substituted in Equation (29) f to estimate the 
angular differential albedo. Primary and secondary energy integrations of 
Equations (24) and (26) are as defined for Equation (20) . 

The shadow shield as a secondary source of both photons and. neutrons., 
ie. scattered and reaction product radiation, has been discussed. Scatter- 
ing to the detector may be accounted for either by the use of the build-up 
factor concept or a Monte Carlo analysis and reaction product radiation 
intensity may -be predicted by the use of Equation (24) or a Monte Carlo 
analysis. In addition to being a second order source, the shield'may also 
be considered as third order source, ie. primary radiation interacting in the 
shield may produce secondary radiation which in turn interacts in the space- 
craft-structure to yield- a third order flux at the detector. 

The detected flux resulting from interactions in the spacecraft structure 
may be predicted in accord -with either Equations (15) , (16) or (26) and 
-discussions thereto, providing the shield -originating flux 0 Q (E a R r^ g^ie^) = 
<£k(E a k,rj s ) , is known and substituted for ^(Eq , r^) . It may be estimated 
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for right-cylindrical shadow shield axially perpendicular emission as, 

?b< E aR - r ls> =f L f<R.ri s . E aR> -S/ E o> .^(E 0 ) n (l-e' XT ) , (27) 

where the function f (R,r ls ,E aR ) , valid for l «r^ s , takes the shield self - 
absorption into account. This function is defined for a cylinder as 


= 2 . fsin X r ls r 

2ir Jo Pj v 


f(R,r ls ,E aR ) 
where 


Y2 © 
Yl 


-n(E a R)(y-yi) 


dy , 


(28) 


yi 

*2 

♦ 


= r ls cos 8 - R cos $ , 

= r ls cos £ + R cos $ , 

o'- “1/ n s Sin R . 

= tt “ Sin (■ 1S R -- ) , 


y 

P 

r ls 


integration variable; a distance , 

integration variable , an angle , 

the shield-to-spacecraft 'area dA‘ distance; 
analogous to r-^ . 


For reasons of clarity the fluxes discussed in this report section are 
schematically summarized in Figure 3 . 


• 2.-3, Code Logic 

During the present report period Code S0SC a spacecraft shield 

optimization study code was designed in the FORTRAN -IV compiler language 
for the NASA-GSFC IBM-3 6 0/91 digital computer. The three major .component 
sub-codes which make up .S0SC: XEST, NUGAM1 and ALB, have been developed , 
debugged and are provisionally operational. Preliminary development of 
codes XEST and ALB was carried out on a CDC-6600 digital computer, accessed 
from the NUS Corporation CDC-Marc II high-speed terminal facility. 
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Preliminary data and resulting observation obtained with the S0SC component 
codes are given in Section 2.4. Code S0SC consists of a main controlling 
program MAIN, a shield thickness prediction program XEST , a com- 
plex geometry scatter flux program ALB, and a Monte Carlo build-up 

factor and albedo calculating program NUGAM1. Development work on 

these programs has focussed on’ gamma photon transport as noted in Section 

1 of this report. Neutron transport has been evaluated according to 'removal' 
(2,13) 

theory 


MAIN executes input and final output operation. Details of the input- 
output are omitted from this report since they are still subject to some change. 
The main programs calls ALB to determine F for Equation (20) . It calls XEST 

to estimate an approximate thickness value, L m j n , for Equation (20) . 

NUGAMT is called to determine build-up factors corresponding to the 
estimated value of L m i ne Using the thus determined values of F and the 
build-up factors, the code solves Equation (18) for the optimum thickness 
L m 4 n by iteration. It will optionally recall NUGAM1 to determine the de- 
viation of the build-up factors for L min from those' for L m i ne . If this 
deviation exceeds a tolerance value the code will reiterate. Code S0SC 
logic is summarized in Figure 4. 


Code ALB determines the angular-energy integrated flux scattered to the 
detector by the spacecraft complex structural components illuminated by 
primary source radiation. Code ALB consists' of an albedo package and 
a generalized geometry package. In its present form the albedo package 
is coded for gamma photons because of the lack of fast neutron differential 
number albedo data. It is proposed that fast neutron data be generated by 
a proposed Monte Carlo code. 

The albedo routines in code ALB determine the scattered energy integrated 
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flux 0 (E ) defined by Equation (13). The main calling program carries 
«£ as 

out the integration over primary source energies . The gamma photon number 
current albedos defined by Equation (15) were obtained from the relationship 


a(E o '0o' 9,< P' t;i ) 


where 


a(E Q ;i) 


g(t) 


= a(E Q , 0,0,0 ,oo;i) .f(0 ) .cos0 .g(t) , (29) 

= a(E;i) . f (0 ).cos0.g(t) , (30) 

o o 

= the angular differential number current albedo 

for gamma photons perpendicularly incident 

(6 =0) and emergent (0=0) from scattering 

material i, ie. 180° (-Q ) backs catter, 

s 

= a function to account for reduced backscattering 
from a material of finite thickness t, 


f(0 ) = a function to account for the albedo behaviour 

o 

with change in 0 

o 

From reference (14) , the azimuthal dependence may be defined in terms of ft , 0 o 
and the total scattering angle 0 S/ as 


-1 


< p = cos 


where 


‘cos0 + COS0 cos 0 
s o 


sin© sin0 


f 


0 = the angle between the incident and emergent photon 

o 

vectors, i. e., the scattering angle. 


The present version of code ALB temporarily assumes 



f(0 o )=Cos0 o , (31) 

g (t) = a(E o ,t;i) /a(E Q ;i) , (32) 

= 1 - e" Ct , (33) 
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where c is a constant such that g(t) = 0 .99 , for t = 2 x (EJ.; X (E Q ) is the 
mean-free-path in material i for photons of energy E 0 . Code ALB uses scatter- 
ing angle 0 S to eliminate the albedo dependence on azimuth <p, in accord with 
the method of reference (14). 

For code preliminary runs experimentally measured values for the perpen- 
dicular number current albedo a (E 0 ;i) were obtained from references ( 7,8). 
For fast neutrons , Monte Carlo data from reference (10) was used for pre- 
liminary evaluations . 

The generalized geometry routine for code ALB requires that the spacecraft 
structure be redefined as simple geometrical shapes, eg. cylinders, rec- 
tangular boxes, prisms, slabs, etc. These shapes allow the flat-sided 
cylindrical spacecraft body, boom-arms, antennae, science platforms, 
etc., to be accounted for. The theory and logic of the geometry code will 
be detailed in the final report to this contract. 


Code NUGAM1 determines the angular energy transport of gamma photons 
in a finite cylindrical shield. It was derived from an existing source 

self -absorption code NUALGAM, developed for NASA-GSFC by NTJS 

Corporation, and described in NUS-536^^. It considers pair-production, 
Compton scattering and photoelectric interaction phenomena . 

In code NUGAM1 , the shield may be composed of a single material of 
simple or complex composition or stratifications ie. "discs" , of simple 
or complex composition, including vacuum. The code is presently designed 
for either an axial point or plane parallel source but may be readily adapted to 
other distributions. Similarly, the code may be readily modified to allow the 
study of annular-cylindrical shields, rectangular slabs, etc. 

.Code NUGAM1 may be used to determine angular differential forward build- 
up and albedos for gamma photons. Preliminary studies with this code, 
which may be either called by code S0SC or used as a separate code, 
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have revealed that the forward build-up factor for the shadow shields is less 
than that to be expected for a large (semi -infinite) shield, in agreement with 
qualitative argument. This is summarily discussed in Section 2.4. 


Code XEST determines the value of L . satisfying equation { 20) and thus 

•min 

obtains W . of equation (17), or the li . of equation (IS). Although 
mm min 

specifically designed for the purpose of shadow shield optimization it is 

coded for larger shields. Code XEST solves equation (20) by the technique 

of iteration. Fbr the first iteration the code assumes a build-up factor of 

unity to determine L . ^ . For the second iteration a Monte Carlo build- 

in m e 

up factor based on L . and calculated by NUGAM1 is used to iterate 

v mm e 

L m in ^ . Iteration, is arrested when 


( L . - L . (h M/L . (h) ^ € , 

mm min |/ mm ' 

e e I e 

w here 

' € = preassigned tolerance , 
h = iteration number . 


(34) 


2.4 Discussion 

Preliminary results obtained with the S0SC component codes are reviewed 
in this report section. They consist of shield thicknesses and weight 
evaluations for assumed typical SNAP-27 RTG source strengths and spectral 
distributions. Example results obtained for the scattering from spacecraft 
structural members are presented ■, and the importance of various factors dis- 
cussed. 


The calculations described in this section are based on a RTG gamma photon 
emission distribution similar to that of the Martin Cronus v ; (thermal load- 
ing of 4100 watts) normalized to a total emission of 1 x 10^ y/ sec. The 
normalization factor has been discussed in Section 1. The RTG fast neutron 
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emission. distribution was based on SNAP-27-1 reported data from- reference 
(1 ); a source emission rate of 5.7 x 10^ n/sec was assumed. For both 
fast- neutrons and gamma photons, axial and axially prependicular emission 
rates were taken as identical, although this is not true in fact. This 
assumption was necessitated by the lack of actual encased RTG source 
data. The assumed source emission spectrum is given in Table I for gamma 
phitons and in Table II for fast neutrons . 

The typical spacecraft for which the calculations were carried out is that 
of Figure 1 , The Outer Planets Explorer, The dimensions of this spacecraft 
were obtained either directly or by scaling, from NASA-GSFC preliminary 
drawings . Figure 5' shows a schematic outline of this spacecraft for the 
discussions in this section. 

Gamma photon cross-section data were taken from references (17) through 
(21) . Neutron cross-section data were obtained from references (13) 
and (22) , 

The Table I source spectrum is that for a 3 year old PuC>2 SNAP-27. Cal- 
culations in reference (16) indicate that the total gamma emission rate will 
increase by a factor of ~ 4 over 18 years. Table III, reproduced from the 
reference, indicates that the energy groups 0.2 to 0.3 , and 2.0 to 3.0 MeV 
are the- most critical, eg. the 2.0 to 3.0 MeV group (2.62 MeV ThC") increases 
by p factor of more than 100 in the first 10 years. This age effect on shield-) 
ing requirements was not studied during this report period because of lack 
of reliable source data . 

The scattering of source gamma photons to the detector by aluminum boom 
tubing proximate to the source (or the detector) was investigated. The 
boom axis was assumed as perpendicularly bisecting a 16“ long unit (ly/sec) 
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source, for E Q = 0.75 MeV. The calculations assumed the boom as coming 
to within 10 cm of the source, although 25 cm is indicated by scaling the 
GSFC drawings. The calculations assumed a line source and a single 
scattering model for a boom tube volume of 0.79 cc/cm. The results of 
the calculations are shown in Figure 6 along with the calculation geometry. 
The detector-incident scattered flux for a single typical boom (boom (1,2) 
in Fig. (5)) is obtained from Figure 6 as 0 Q (E^) = 1.3 x 10 ^ y/cm^-s ec. 
Assuming boom (1 ,2) in Figure 5 as typical, where r^ ranges from 25 to 
80 cm, ie. a length of 55 cm, the total scattered flux from 6 such booms 
(6 per side) and 4 RTG sources would be ~ 3 x 10 y/ cm z -sec. 

The use of the single scattering model as opposed to the albedo method for 
boom structure calculations, was investigated. A 1 cm length of aluminum 
boom was considered as located 100 cm from both a' point source and a 
detector, with source and detector 141 cm apart. For E 0 = 0.75 , a unit 
source (1 7/sec) and a boom volume 0.79 cc/cm, the calculated fluxes 
were 

s (Eq) = 3.6 x 10 12 , y/cm* 2 -sec, 

0 ass (Ej = 6.4x10 12 , y/ cm 2 -sec. 

The difference -in these calculations diminishes if it is assumed in the 
case of the albedo result that photons penetrating the relatively thin (0.04") 
tube frontal wall may be backscattered from the tube interior wall surface 
'ie, if the wall thickness is doubled for the calculation, then 

-12 2 

&as a? = 7,0 x 10 / y/cm -sec. 

This last result is only slightly greater than $ ass (E^) . The conclusion may 
be drawn that the single scattering result is a reasonable and valid approxi- • 
mation. 
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A. similar calculation for a large cylinder, in the physical position of the 
spacecraft cupula, was carried out for iron as the material, dimensions 
20" x 38.5" x 0.5" (dia . x It. x thickness) and incident energy E Q = 0.75 
MeV . The calculated fluxes were 

<p as (E^) = 2.4 x 10 9 y/cm 2 -sec, 

<p a ss (Eg) = 2.7 X 10 9 y/cm 2 -sec. 

■For aluminum and. a thickness of 1 inch, the singley scattered flux was 

^ass = 1.8 x 10 9 y/cm 2 -sec. 

Taking an actual cupula 'flat' , A, in Fig. (5) , with dimensions 13.2" x 25.5" 
(width x It.) and an assumed aluminum thickness of 1" , the albedo technique 
flux was determined as 

0 as (EJ = 1.8x10 10 , y/cm. 2 -sec. 

Since (p^ Q (E^ ~'0 O!SS (E ) / the scattered and detector incident flux from 
the twelve cupola 'flats’ would be ~ 2 x 10 9 , which is in good agreement 
with the large cylinder approximations; this assumes no attenuation of the 
'flat' -transmitted radiation. 

The above calculations were based on a single RTG source of 1 photon/sec, 
and E Q = 0.75. Taking the source as 10® y/sec, the detected scattered flux 
would approximately equal the contract criterion, C ^10 particles/cm 2 -sec. 
However-, recent calculations using the photon spectrum of Table I indicate 

(E ) 0.1 , 7 /cm 2 -sec, 

a 

for the large iron cylinder model , and 

2 

ch (E ) =“ 0.15, y/cm -sec, 

a 
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for the flat aluminum surface, . Although these values are less than the 
single energy calculated values, they indicate that spacecraft structure 
scattering may in some problems be significant relative to the criterion C. It is 
proposed to study this problem further in the remainder of the contract period . If 
this photon scattering is considered significant enough to require shielding 
then such shielding would be a minimum weight if designed to attenuate the 
scattered photons as opposed to the primaries. The scattered photon energies 
will generally be in the range, 0.15 to 0 f 5 MeV, and thus readily attenuated 
by a thin shield located at the detector. 

The foregoing calculations were obtained with the S0SC component code , 

ALB. Figures (7 ) and (8 ) are sodium-iodide scintillation detector spectra 
reproduced from reference ( 7 ) which indicate the ratio of single to multiple 
scattering as a function of t and Q q , respectively, for graphite as the back- 
scattering material. The prominent peak in these spectra is the single scatter 
distribution . 

Table IV gives example optimum results for shield thicknesses and weights 
as predicted' by code XEST, for lithium -h ydride , aluminum, iron and lead. 

The calculations assumed the source spectra of Tables I and II. 

Build-up factors of 1.05 , 1.3 and 1.5 were used for both gamma photons 
and neutrons in these preliminary calculations. The choice of these values 
is discussed below. The calculations do not include reaction product 
radiation. 

The preliminary photon number build-up factors used in XEST were obtained 
from the Monte Carlo code NUGAM1. A point source axially located 8 cm 
from the shield was used throughout. Investigations are now being carried 
out for radiation axially- parallel incident. Early determinations were based 
on jljL ~ 2 .0 from which build-up factors of ~ 1 . 05 resulted. Recent deter- 
minations have indicated that higher values are likely since, as Table IV 
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indicates, juL~4. For E q = 1.25,, the NTJGAM1 point source code version predicts 
build-up factors ~1 . 3 , for juL~4 . 0 . 

Examples of energy-integrated scatter photon angular distributions are presented 
in Figure 9 . The angular categorization is referenced to the shield geometric 
center. The distributions were determined by the Monte Carlo subprogram: 
NUGAM1. The current code version determines the photon distribution to which 
the detector is specifically exposed, i.e. , photons escaping from the shield 
sides are considered' as not intercepting the detector. 

Table V gives sample number build-up factors as determined by code NUGAM1 

for the noted materials and photon energies. They are based on 10,000 source 

photons uniformly distributed in the source -shield solid angle subtended by the 

circular shield face. The build-up factors in this Table must be considered 

as high since they include some shield side escape i.e. , the forward cate- 

o 

gorizing solid -angle had a half -angle of 10 . 

Further progress with respect to the foregoing will be reported in the monthly 
progress reports and the contract final report. 
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3. WORK PROGRAM FOR SECOND TERM 


Throughout the second and remaining term of this contract, the work program 
will be carried to completion. Code S0SC will be integrated, debugged and 
made operational on the NASA/GSFC IBM-360/91 digital computer. Code 
S0SC will be run to- determine the optimum shield for the spacecraft example 
of Figure 1. The preparation of a-detailed user’s manual and code description 
will complete the work scope outlined under Tasks I and II of this contract. 

The design of a suitable and appropriate laboratory experiment will be carried 
out to verify the code S0SC shielding predictions. This result of this design 
will be incorporated into the code user/description manual, which will repre- 
sent the contract final report. 
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FIGURE 5 

SCHEMATIC DRAWING OF SPACECRAFT 
SHOWING SPECIFIC STRUCTURAL MEMBERS 
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FIGURE 6 


BOOM SCATTER INTENSITY AS A FUNCTION 
OF SOURCE-TO-BOOM DISTANCE 




400 Min; X = Material Thickness 



FIGURE 7 

Nal(Tl) Scintillation Spectra Showing 
Gamma Photon Backscatter Distribution as a 
Function of Scattering Material Thickness 



Cjs 137 ; Graphite 

400 Min; 9 = 80-(9 +9) 
• s o 



8 (a) Perpendicular Incidence: 

Spectral Distribution as a Function 
of Scatter Angle 9g for X 0 = 0.668 MeV 

and 9 =0° 


Cs^ 37 ; Graphite 

400 Min; 9 =180 -(9 +6) 
s o 



8 (b) Slant Incidence: 

Spectral Distribution es a Function 

of Scatter Mel® 9 S f or E c = 0.662 MeV 

and 9 =40° 
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FIGURE 8 


Nal(Tl) Scintillation Spectra Showing 
Gamma Photon Scatter Distribution as a 
Function of Incident -Angle , 0 Q 
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FIGURE 9 


ENERGY-INTEGRATED ANGULAR DISTRIBUTION 
OF SHIELD SCATTERED PHOTONS 






TABLES I to V 
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TABLE I 


Assumed SNAP-27 Gamma Photon Emission Spectrum 

( 16 ) 

Based On Martin Cronus Data' 7 


Energy . ' 

Interval 
(MeV) ' 

Assumed 

Energy 

(MeV) 

Photon Emission 
Rate 

(v/MeV-sec) 

0.044 

- 0.2 

0.15 

6.54 

X 

10 7 

0.2 - 

0.3 

0.24 

6.87 

X 

10 7 

0.3 - 

0.4 

0.311 

7.88 

X 

10 7 

0.4 - 

0.5 

0.414 

7.86 

X 

10 7 

0.5 - 

0.6 

0.5 83 

8.03 

X 

10 7 

0.6 - 

0.7 

0.650 

7.33 

X 

10 7 

0.7 - 

0.8 

0.766 

2.26 

X 

10 8 

0.8 - 

0.9 

0.851 

1.40 

X 

10 8 

0.9 - 

1.0 

1.00 

6.65 

X 

10 6 

1.0 - 

1.2 

1.10 

1.87 

X 

10 7 

1.2 - 

1.4 

1.40 

1.73 

X 

10 7 

1.4 - 

1.6 ‘ 

1.59 

7.71 

X 

10 6 

1.6 - 

1.8 

1.63 

1.63 

X 

10 7 

1.8 - 

2.0 

1.90 

1.54 

X 

10 7 

2.0- - 

3.0 

2.61 

1.07 

X 

10 8 

3.0 - 

4.0 

3.50 

3.03 

X 

10 s 

4.0 - 

5.0 

4.50 

9.10 

X 

10 3 

5.0 - 

6.0 

5.50 

3.01 

X 

10 4 
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TABLE II 


SNAP-27 -1 Fast Neutron Emission Spectrum 


d) 


Energy 

Interval 

(MeV) 


Neutron Emission Rate 
(n/MeV-sec) 


0 - 1.0 

1.0 - 3.0 

3.0 - 4.0 

4.0 - 5.0 

5.0 - 6.0 

6.0 - 8.0 
8.0 - 10.0 


5.06 x 10 7 
4.16 x 10 7 
1.79 x 10 7 
3.56 x 10 6 
1.90 x 10 6 
4.50 x 10 5 
1.14 x 10 5 


Energy Integrated Emission Rate = 5.7 x 10 


7 


n/ sec. 
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TABLE III 


ReDroduced from Reference (16) 

2 

Gamma spectra, photons /cm -sec for RTG at various times 
after plutonium separation (Normalized to 1.0 for 
energy interval 0. 5-1.0 MeV at 18 years) 


Energy 

(MeV) 

0 year 

1 year 

5 year 

10 year 

18 year 

0. 04 

- 0. 5 

0. 12 

0. 16 

0. 56 

0. 90 

1. o 

0.5 

- 1. 0 

0.23 

0. 25 

0. 30 

0. 33 

0. 35 

1. 0 

- 2 

0.096 

0. 096 

0. 096 

0. 096 

0.096 

2 

- 3 

0.003 

0. 018 

0. 18 

0.32 

0. 37 

3 

- 5 

0. 0003 

0. 0003 

0. 0003 

0. 0003 

0. 0003 

5 

- 7 

0. 00003 

0. 00003 

0. 00003 

0. 00003 

0. 00003 
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TABLE IV 


Optimized Shield Thicknesses And Weights 
Based On Provisional Build-Up Factor Data* 




B = 

1.05 

B = 

1.3 

B = 

1.5 

Shield 

Material 

. density, p 
(g/cc) 

length , L 
(cm) 

weight, W 
(kg) 

length, L 
(cm) 

weight, W 
■ (kg) 

length , L 
(cm) 

weight , W 
(kg) 

lithium - 
hydride 
(Li H) 

0.82 

56.28 

1.46 

60.90 

1.58 

64.1 

1.66 

aluminum 

(Al) 

2.699 

18.3 

1.57 

19.8 

1.69 

20.7 

1.77 

iron 

(Fe) 

7.87 

7.44 

1.85 

8.10 

2.02 

8.54 

2.13 

lead 

(Pb) 

11.35 

5.80 

2.08 ■ 

6.54 ■ 

2.35 

7.08 

2.54 


*Right -Cylindrical Shield of Radius 1.25 inches. 



TABLE V 


Provisional Build-Up Factors Predicted By Code NUGAM1 


Shield 

Material 

length, E 
(cm) 

Photon 
Energy, E Q 
(MeV) 

Number of 
Mean-Free 
Paths, ji(E 0 )L 

Flux at Dete.c.tor 
(y/ cm^) 

Unscattered Scattered 

. Build-Up 
Factor , 

B 

LiH 

29.1 

‘ 1.25 

1.37 

.0450 

.000615 

1.025 

LiH 

29.1' 

0.25 

2.7 

.0115 

.00121 

1.05 

A1 

67.7 

2.5 

7.05 

.000153 

.000868. 

6.67 

A1 

27.5 

1.25 

4.07 

.00301 

.00102 

1.34 

A1 

1.4. 8 

0.75 

2.81 

.00112 

.0112' 

1.10 

Pe 

4.89 

0.75 

2.64 

.001,24 

.0126 

.1.10 

Pb ' 

5.28 

1.25 

3.46 

.00104 

; 00554 

1.188 
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APPENDIX 

INTERACTION PHYSICS REVIEW 


I GAMMA PHOTON INTERACTION PHENOMENA 


In their passage through a medium, photons interact with the electrons 

and'nuclei of atoms in their path. These phenomena form the basis both 

for their detection and for the deposition of their energy. A discussion 

of the kind and effect of these- interactions as they pertain to the present 

(23-27) 

work is given in this report section ' . 


There are four kinds of basic gamma photon interaction processes 
of which only two are relevant in the present work, namely: 


(27) 

i 


a) interaction with atomic electrons , 

b) interaction with the electric field 
surrounding nuclei or electrons . 

The effect of (a) may be either scattering or absorption; the latter is 
the Photoelectric Effect. The scattering may be either one of the two 
types: 


1 . Compton inelastic scattering (incoherent) , or 
2 . Rayleigh elastic scattering (coherent) . 


The effect of (b) is the disappearance of the photon and the creation 
of an electron- pair; this phenomenon is referred to as the Pair Production 
Effect. 

A brief discussion of these four microscopic phenomena, and their 
macroscopic attenuating effect on a beam of photons, follows under 
the headings: 


A) Photoelectric Effect 

B) Compton Scattering 

C) Rayleigh Scattering 

D) Pair Production 

E) Attenuation 
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A. Photoelectric Effect 


At relatively low photon energies the most probable effect of an interaction 

is absorption of the incident photon by an electron of the traversed medium 

followed by ejection of that electron and emission of either characteristic 

X-rays or “Auger electrons" as explained below. This phenomenon, called 

the- Photoelectric Effect, results in the complete disappearance of the inci- 
( 23 ) 

dent photon' 

In order that total absorption may take place, and momentum be conserved, 
the interacting electron must be initially bound, in which case . the residual 
atom recoils. The most tightly- bound electron, with respect to the incident 
photon energy, has the greatest probability of absorbing the photon. The 
interaction cross-section is a maximum when the photon energy E y, is just 
equal to the electron binding energy; it decreases gradually as Ey increases, 
and decreases sharply as Ey decreases. The most tightly- bound electron in 
an atom is in the K- shell; it accounts for in excess of 80 % of the photo- 
electric absorptions, with the L-shell accounting for most of the remainder. 

The energy of the ejected electron, or photoelectron as it is usually called, 
is given by: 



E o = hu - E (MeV) 

e o eb 

where 

- hu = incident photon energy (MeV). 
o. 

• E , = electron binding energy (MeV) . 
eb 


( 1 ) 


The energy E . is carried away from the atom by radiation emitted as the 
eo 

inner shell vacancy is filled by an outer shell electron, such radiation is 
referred to as Characteristic X-rays . If the X-rays interact with an outer 
shell electron as they leave the atom, they will be absorbed and the 
absorbing electron emitted instead — an Auger electron. The nuclear 
decay processes of internal conversion and electron capture may also 
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lead to the emission of characteristic X-rays; so also will the absorption 
of beta particles . 


The photoelectric effect does not lend itself easily to explicit theoretical 

calculation. Determinations of its cross-sections are usually based on a 

combination of empirical treatments which vary according to the energy 

range under consideration. It is the practice of most researchers to make 

(17-21) 

use of tabulations for or , the photoelectric cross- section' 

1 iL 


B. Compton Scattering Effect 

As the wavelength of gamma photons decrease and their photon energy 

increases , their behaviour tends towards that of a particle and their 

identity with a wave diminishes . The region of this transition corresponds 

to the Compton scattering "threshold." This threshold, not sharply- 

2 

defined, is entered upon gradually as h u m c , (=0.51 MeV), where 

o o 

m equals the rest mass of the electron. Viewed as solid bodies, the 
o 

photon and the electron have comparable "masses." As the Compton effect 

(23) 

becomes significant, the photoelectric effect significance diminishes' 

Compton scattering may be considered as an inelastic collision between an 

incident photon and a "free" electron of the medium; the collision is 

analogous to that of billiard ball mechanics . The electron may be thought 

of as free to recoil on the basis of hn » E , , as a result of which the 

o eb 

incident photon may, transfer a portion of its momentum and energy. The 
consequence of the collision is a scattered photon of energy hu^, travelling 
in a new direction and at an angle 0 with the original photon direction , and 
a recoiling electron of -energy E g making an angle $with the incident photon 
direction. 0 max = 180°; |tmax = 90°. 

The angular and energy relationships of these statements may be expressed 
as; 
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MeV. 


( 2 ) 



i + — (1 - cos e) 


m c 
o 



.E = hu , MeV. 

y 

/ * \ 0 

cot; 1) = (1 + a ) tan -7- 
o l 


(3) 

(4) 

(5) 


where 

E = incident photon energy, MeV. 

^ o 

E = scattered photon energy, MeV. 
y l 

0 = angle between incident and scattered 

photon directions . 

ip - angle between incident photon and 
recoil electron directions . 


For convenience in the -remainder of this section the following conventional 
short form is used: 




me me 

o o 


(6) 


The differential "collision" cross-section for the scattering of photons 
into a given solid angle dQ at a particular angle 0 is given by the -Klein- 
Nishina formula as 
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dcr = 


dQ 


( 1 + cos 0) 




jj + a' (1 - cos 0)J 


2 

a (l-cos0) 2 


x <l+- 


J 


(1 +cos 0 ) {j. +a' ( 1 - cos 0)J 


where 


d ct 
dQ 


differential cross-section, cm /electron 

2 7T Sin 0 d 0, the differential solid angle 

2 

"classical electron radius , — — - 


m c 
o 


(7) 


2.818 x 10 13 cm. 


Equation (7) assumes the incident photons .to be unpolarized. It 
indicates that for large a ’ q , scattering is predominantly in the forward 
cone. As a ■ 0, and cos 0 — 1 , we see 


do 


~ — ( 1 + cos 2 0 ) da 


(8) 


From Equation (2) equation (7) may be rewritten in terms of energy 


d q- 
da' 


irr 


i- + JL + i 


a 


a a 
o 


i 2 i 2 ii 

a a , a a 


o 


o 1 


a. a , 

+ -°-+-k 

l t 

a . a 


0) 


for 


a 


a a . 
o 1 


(1 + 2a ) 
o 
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The integration of Equation (9} over all scattered energies yields 

the total Compton scattering cross-section per electron, a , : 

cs 


cs 


= 2irr 


9 ! 1 + a ' 

2 ] o 


O' 


o 


2a (1 +a ) 

- — - ln(l + 2a' ) 

1 + 2a’ ° 


+ 


+ 


In (1 + 2o: ) 


1 + 3tt 


2a 


o 


(1 + 2 a Y 

o 


( 10 ) 


The total Compton scattering cross-section per atom is given by Z • a , 

cs 

where Z is atomic number. 


C . Rayleigh Scattering Effect 

In Compton scattering the atomic electrons are assumed to be unbound. 
This assumption is only valid at photon energies which are large with 
respect to the electron binding energy. A low energy photon may be 
elastically scattered by a tightly bound atomic electron, with the atom 
as a whole absorbing the recoil momentum. A bound electron has a 
"mass" which is equivalent to that of its atom. The energy transferred 
to the atom is small, and so the scattered photon proceeds with a 
relatively unaltered energy and only a slightly altered direction. This 
effect is known as the Rayleigh or small-angle scattering effect. 

Since all the electrons in a given atom behave similarly, Rayleigh' 

scattering is coherent. Because all the atoms of a given solid may be 

packed regularly, the effect may extend to the electrons of different 

atoms. When the scattering angle, 0 =** 0, the scattering will be in 

K 

phase, i.e. constructive interference. As 0 increases the tendency 

K 
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The integration of Equation (9) over all scattered energies yields 

the total Compton scattering cross-section per electron, a , : 

cs 


a = 2irr 
cs ( 


2 I 1+a o 


a 


. 3 


2 a ( 1 +.a } 

— _ in (1 + 2a' ) 

1 + 2a' ° 


+ 


In (l + 2 a ) 
o_ 

2a 


1 + 3a 


(1 + 2a' Y 
o 
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The total Compton scattering cross-section per atom is given by Z * cr , 

cs 

where Z is atomic number. 


C . Rayleigh Scattering Effect 

In. Compton scattering the atomic electrons are assumed to be unbound. 
This assumption is only valid at photon energies which are large with 
respect to the electron binding energy. A low energy photon may be 
elastically scattered by a tightly bound atomic electron, with the atom 
as a whole absorbing the recoil momentum. A bound electron has a 
"mass" which is equivalent to that of its atom. The energy transferred 
to the atom is small, and so the scattered photon proceeds with a 
relatively unaltered energy and only a slightly altered direction. This 
effect is known as the Rayleigh or small-angle scattering effect. 

Since all the electrons in a given atom behave similarly, Rayleigh' 

scattering is coherent. Because all the atoms of a given solid may be 

packed regularly, the effect may extend to the electrons of different 

atoms. When the scattering angle, 0^ 0, the scattering will be in 

phase, i.e. constructive interference , As 0 increases the tendency 

K 
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is towards destructive interference and so the scattered photons will 
be found concentrated mainly in a narrow forward cone, and to a lesser 
extent in other discrete directions. This may be realized from consid- 
eration of the photon wavelength, and the atomic radius, analogous to 
(24) 

Bragg reflection . This behaviour differs from Compton scattering, 
where the independence of the electrons precludes the liklihood of 
interference. 

The transition from Rayleigh scattering to Compton scattering is smooth 

with increasing energy, E r . The Rayleigh scattered photon does not 

1 o 

have a unique energy as a function of scattering angle, having instead 
an energy distribution peaked at a value close to that given by Equation (2) . 


D. Pair Production Effect 

At photon energies of approximately 1 . 0 MeV the predominant interaction 

phenomenon is Compton scattering. As E., is increased considerably 

to 

above this energy the photon may interact with the electric field surrounding 
either a nucleus or an electron. The photon will be absorbed and replaced 
by a pair of electrons, a positron and a negative electron. This effect is 
called Pair Production. 

The cross-section for pair production in the field of an orbital electron 

2 

is negligible until E v $-4 me , (= 2.04 MeV). Nuclear pair production , 

' o o 

however, has a cross-section which begins at the photon threshold energy 
2 

2m c , (= 1.02 MeV), and increases rapidly thereafter, 
o 

The electron pair share and carry away the energy in excess of that 
required for their creation, as kinetic energy; this may be expressed as 
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( 11 ) 


( E _ + E + ) = E - 2 m c 
e y o 

o 

The free positron is quickly annihilated by a negative electron after its 

kinetic energy has been dissipated. The annihilation yields a randomly 

2 

oriented pair of back-to-back photons , each with an- energy of m c . 

o 

The electron pair are distributed mainly in the forward direction with the 

2 

average angle of "deflection" being expressed by m c /E . 

2 2 2 o e 

For 2 m c < E < 4m c . a orZ . 

o y o pp 

o 


E. Attenuation 

The passage of a beam of photons through a medium is characterized by 
their interactions with the atoms of that medium. This leads to a reduc- 
tion in the number of uncollided primary photons at a depth . The reduc- 

(2 23 24) 

tion is referred to as the attenuation of the incident photon number* ' ' 

The discussion on interactions has shown that the total microscopic 

energy dependent cross-section for a particular interaction process 

occurring, is given by either or , Z.cr , ct do or cr . The total 

cross-section a E v ), for "some" process occurring, is then-given 
TOT / o 

by the sum of the partial cross-sections as 

CT TOT (E y > = a TOT = g PE +Z ' a CS +a RS + °PP ' <=m 2 /atom) (12) 

o 

from which a total macroscopic cross-section per cm of path may be 
defined as 

M TOT (E y o } = ^TOT =Na TOT' (cm ” 1) 
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The free positron is quickly annihilated by a negative electron after its 

kinetic energy has been dissipated. The annihilation yields a randomly 

2 

oriented pair of back-to-back photons, each with an energy of m^c • 

The electron pair are distributed mainly in the forward direction with the 

2 , 

average angle of "deflection" being expressed by m c /E . 

o e 

2 • 2 2 
For 2 m c < E < 4m c , 0 - ccZ . 
o v o pp 

o ^ 


E. Attenuation 

The passage of a beam of photons through a medium is characterized by 
their interactions with the atoms of that medium. This leads to a reduc- 
tion in the number of uncollided primary photons at a depth. The reduc- 

(2 '23 24) 

tion is referred to as the attenuation of the incident photon number ' ' 

The discussion on interactions has shown that the total microscopic 

energy dependent cross-section for a particular interaction process 

occurring, is given by either cr __ , Z.cr , a n _ or cr . The total 

PE Ob Kb r f 

cross-section o ( Ey ) , for "some" process occurring, is then.given 
by the sum of the partial cross-sections as 


a TOT (E y 1 = ff TOT = “pE + Z ‘°cs + °RS + °PP ' (cm M ° m) 


( 12 ) 


from which a total macroscopic cross-section per cm of path may be 
defined as 

^TOT ^ E y ^ = ^TOT~ WCT TOT' ^ 

o 


(i3) 


- 53 - 



and similarly 


“PE = N °PE : ^CS = Nz ff CS ■' "RS = N ct RS ; »Vp = N °PP <M> 

where 

p x A 

N = — ; . (atoms present/cc) ( 15 ) 

assuming one type of, atom only. 

p = density of medium, (gm/cc). 

-23 

A^— Avogadro number, (atoms/mole); 6.023 x 10 
M = atomic or molecular weight. 


The inverse of Equation is defined as the mean free path, A , 

TOT 

for a photon prior to interaction , i . e . 


X 


TOT 


, (cm) 


TOT 


(16) 


Similarly X , X , X , and X may be defined from Equation (14) 
Kiii ob Ku ii 

The total macroscopic cross-section is generally referred to as the total 
linear attenuation coefficient. 


The number of normally incident photons per cm -sec in a parallel beam 
which penetrate a thickness x of a homogeneous medium without inter- 
action, is given by the exponential law as 


$ (x) = $ (0)e X x (-y/cm^-sec) (17) 

where 

2 

(0) = number of photons/cm -sec incident at x=0. 

2 

$ (x) = number of photons/cm -sec emerging at x. 

= cross-section appropriate to the interaction 
effect under consideration; energy dependent. 
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Thus, it follows that the probability of uncollidqd photon transmission 
through a thickness x is given by p , where 


P = 


0 (x) 

0 ( 0 ) 


^TOT X 


= e (18) 

Equation (18) is in agreement with Equation (4) of Report-Section (2.2) 
Similarly the probability of some kind of interaction occurring in 
path length x is given by 


^TOT X 

€ tqt (x) = (1-p) = (1 - e iUi ) 


(19) 


Equations , and assume a normally incident 

parallel or collimated, photon beam. In practice this can only be nearly 

achieved by narrow- geometry restrictions . For the case of a poorly 

collimated or uncollimated photon beam it is necessary to introduce a 

factor to express the increase of the photon number flux at x over the 

valued predicted by Equation (3.4.17). This factor, known as- a number 

(2 3 241 

Build- Up Factor, B(x), may be defined as x ' ' ' 


B (x) = 


Total Number Flux at x 


£ 1.0 


Uncollided Flux at x 
which is in agreement with Equation (5) of Report-Section (2.2) 


( 20 ) 


Energy and dose build-up factors may be similarly defined. 


That /j, is a function of both incident photon energy, E y , and the properties 

of the traversed medium is apparent from the discussions of sub sections 

(A-D) and equations (12) to (.15) . It follows then, that any property, 

including B (x) which is dependent on \i, is similarly dependent on E v , 

' o 
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II FAST NEUTRON INTERACTION PHENOMENA 

Neutrons interact with the nuclei of traversed matter through the mechanism 
of nuclear force. An interaction is generally referred to as either having 
scattered or absorbed the incident neutron. The probability of either scattering 
or absorption varies as a function of the incident photon energy and the 
atomic number of the target nuclide , the dependence on atomic number being 

general and such that for each isotope there is a unique probability, or 

.. ( 2 , 28 , 29 ) 

cross-section 

Perhaps the most convenient and systematic manner of describing neutron 
interactions consists of invoking the "compound nucleus" concept. According 
to this concept all interaction modes result in the formation of an intermediate 
reaction product - a compound nucleus - formed by an absorption of the incident 
neutron. Symbolically this may be represented as 

i n ° + / ><z^ +1 >* 

o 

where, ]_n is the incident neutron, the target nucleus before interaction 
A + 1 ^ 

and x the compound nucleus formed by the interaction. The asterisk 

La 

denotes that the compound nucleus will, in general, be left in an excited state 
for a finite period of time. The energy of the compound' nucleus includes 
both the binding energy and/or part of the kinetic energy of the incident 

fa 

neutron. This energy excess over that of Z X , distributed in a complex 

fashion among the nucleons, will cause the compound nucleus system to seek 

its state of lowest permissible energy in a' characteristic "relaxation time" , 

-20 -12 

typically about 10 to 10 seconds. 

The laws of quantum mechanics allow only those reactions to take place 
which obey certain rigid energy and momentum relationships . This obedience 
is observed with respect to the available excitation energy of the 
incident neutron and the distribution of energy levels in the target nucleus , 
as well as with respect to the symmetry requirements of the interaction, 
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and Jx the compound nucleus formed by the interaction. The asterisk 
denotes that the. compound nucleus will, in general, be left in an excited state 
for a finite period of time. The energy of the compound nucleus includes 
both the binding energy and/or part of the kinetic energy of the incident 

A 

neutron. This energy excess over that of , distributed in a complex 

fashion among the nucleons, will cause the compound nucleus system to seek 

its state of lowest permissible energy in a characteristic "relaxation time" , 

-20 -12 

typically about 10 to 10 seconds. 

The laws of quantum mechanics allow, only those reactions to take place 
which obey certain rigid energy and momentum relationships . This obedience 
is observed with respect to the available excitation energy of the 
incident neutron and the distribution of energy levels in the target nucleus , 
as well as with respect to the symmetry requirements of the interaction, 


55 - 



e. g. parity, baryon number , charge conservation, statistics. 


In accord with the compound nucleus concept, scattering and absorption of 

an incident neutron may be defined in terms of whether or not the compound 

nucleus emits a neutron during de-excitation. Further and more important, 

scattering may be separated into two kinds - elastic and inelastic. In 

elastic scattering the compound nucleus emits a kinetic energy degraded 

-20 

neutron in a very short relaxation time <10 seconds, and is itself left in 
exactly the same internal energy state as before the interaction. In inelastic 
scattering the compound nucleus emits a neutron of partially or totally degraded 
kinetic energy and is itself left in an internal energy state above that of 
X A ; the excess energy is evolved by emission of one or more gamma photons. 

2j 

The degradation of the neutron kinetic energy by scattering is referred to 
as thermalization . The compound nucleus may be de-excited by emission of 
particles other than neutrons, such as alphas and betas accompanied by gamma 
photons in which case neutron absorption is said to- have resulted. 

The interaction processes reviewed above are summarized: 


(i) Elastic Scattering, (n,n) A neutron of reduced kinetic energy 

is emitted by the short-lived compound nucleus which is left in an unexcited 
state. Kinetic energy is transferred to recoil the target nucleus. In the 
case of hydrogen target nuclei, for which energy transfer is a maximum, 
recoil protons result. 


(ii) Inelastic Scattering, (n,n‘) A neutron of reduced kinetic 

energy is emitted by the compound nucleus. The compound nucleus is de- 
excited by gamma photon emission. 

(iii) Radiative Capture, (n,y) The compound nucleus formed by 

absorption of an incident neutron is de-excited by relatively high energy 
gamma photon emission. 
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(iv) Charged Particle Emission', (n,p) , (n,d) , (n,ry) The compound 

nucleus formed by absorption of an incident neutron is de-excited by emission 
of a charged particle such as a proton, deuteron or alpha, accompanied by 
gamma photons. 


(v) Fission, (n,f) The compound nucleus formed by absorption of 

an incident neutron breaks into two ionizing fission fragments, and one or more 

# 

energetic neutrons , .accompanied by gamma photon emission. Fission is most 
probable in heavy nuclei of odd mass number and less so in heavy nuclei . 
of even mass number. It can occur either as the result of an externally 
incident neutron or as a consequence of quantum mechanical leakage through 
the Coulomb barrier, ie. , spontaneous fission. 


(vi) Other Reactions , (a>,n) , (v,n) Two interaction processes which 

'give rise to neutron emission and thus which must be identified are (a,n) and 

(•y,n) phenomena in the plutonium -oxide source and its immediate environment. 

The first reaction proceeds when the energy of an alpha particle exceeds the 

energetic threshold and Coulomb repulsion barrier for the reaction, and is 

thus significant only for plutonium alphas (~5.5 meV) incident on light-target 

nuclei such .as beryllium. The second reaction, photo-neutron formation, 

results from the interaction of high energy gamma photons with light nuclei 

such as are present in plutonium -oxide as impurities. High energy photons are 

2 08 

present in the PuO source through (n ,f) reactions and the decay of the 'll 
2m3 6 

daughter of the Pu isotope present in plutonium oxide . 
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